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I. INTRODUCTION 

Quantum Chromodynamics (QCD) is widely believed 
to be the theory of the strong interactions. Investiga- 
tions on its rich vacuum structure and how the QCD vac- 
uum can be modified in extreme environment are among 
the major theoretical challenges in modern physics. It 
is in particular an interesting topic to study how non- 
perturbative features of QCD are affected by thermal ex- 
citations at high temperature T and/or by baryon-rich 
constituents at large baryon (quark) chemical potential 
fiq. Such research on hot and dense QCD is important 
not only from the theoretical point of view but also for 
numerous applications to the physics problems of the 
Quark-Gluon Plasma (copiously produced in relativistic 
heavy-ion collisions) , ultra-dense and cold nuclear/quark 
matter as could exist in the interior of compact stellar 
objects, and so on. 

The most intriguing non-perturbative aspects of the 
QCD vacuum at low energy are color confinement and 
spontaneous breakdown of chiral symmetry. In recent 
years our knowledge on (some parts of) the QCD phase 
diagram has increased noticeably because of significant 
developments of the lattice QCD simulations (see 
01 for several examples and see also references therein). 
At zero fiq, except for some reports 0, the numeri- 
cal simulations have almost established that two QCD 
phase transitions (crossovers) take place simultaneously 
at nearly the same temperature; one for quark deconfine- 
ment and another for restoration of chiral symmetry (the 
latter being always broken because of finite bare quark 
masses, strictly speaking). It is still however under de- 
bate whether two crossovers should occur at exactly the 
same temperature, however. 

Once a finite is turned on, the Monte-Carlo simu- 
lation in three-color QCD on the lattice cannot be per- 
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formed straightforwardly because of the (in) famous sign 
problem To overcome this problem several tech- 

niques have been developed such as the multi-parameter 
reweighting method @, Taylor expansion Q, density of 
state method Q , analytical continuation from the imag- 
inary chemical potential the complex Langevin dy- 
namics [1^, etc. 

In addition to hot and dense QCD with T and fiq, 
the effect of a strong magnetic field B on the QCD vac- 
uum structure is also a very interesting subject. It would 
be of academic interest to speculate modification of the 
vacuum structure of a non-Abelian quantum field theory 
under strong external fields. Besides, more importantly, 
this kind of investigation has realistic relevance to phe- 
nomenology in relativistic heavy-ion collisions in which 
a strong magnetic field is produced in non-central col- 
lisions [III, In particular, the results obtained by 
the UrQMD model [ijl show that eB created in non- 
central Au-Au collisions can be as large as eB « 2m^ 
(i.e. B ^ lO^^Gauss) for the top collision energy at RHIC, 
namely ^/s^ = 200 GeV. Moreover, an estimate with 
the energy reachable at LHC, ^/s^ « 4.5 TeV, gives 
eB w 15to^ for the Pb-Pb collision according to Ref. [l^ . 
Hence, heavy-ion collisions provide us with a most in- 
triguing laboratory available on the Earth in order to 
study the effect of extremely strong magnetic fields on 
the QCD vacuum. 

Concerning the (electromagnetic) magnetic field effect 
on the QCD vacuum structure, there have been many 
investigations and it has been recognized that B plays a 
role as a catalyzer of dynamical chiral symmetry break- 
ing [l3l-[l6|. The QCD vacuum properties have been also 
studied by means of so-called holographic QCD mod- 
els The relation between the dynamics of QCD in a 
strong magnetic field and non-commutative field theories 
is investigated in Ref. [13] . 

A phenomenologically interesting consequence from 
the strong B in heavy-ion collisions is what is termed the 
Chiral Magnetic Effect (CME) The underlying 

physics of the CME is the axial anomaly and topological 
objects in QCD. Analytical and numerical investigations 
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have demonstrated that the sphaleron transition occurs 
at a copious rate at high temperature unhke instantons 
that are thermaUy suppressed [lO, [2l|. Sphalerons are 
finite-energy solutions of the Minkowskian equations of 
motion in the pure gauge sector and they appear not only 
in the electroweak theory but also in QCD They 
carry a finite winding number Qw which is defined as 

= y d^x Tr[FF] , (1) 

where F and F denote respectively the field strength ten- 
sor and its dual. Sphalerons connect two distinct classical 
vacua of the theory with different Chern-Simons numbers 
in Minkowskian space-time. It is possible through the 
coupling with fermions in the theory to relate the change 
of chirality, Nn — , to the winding number by virtue 
of the Adler-Bell- Jackiw anomaly relation, 

{Nr - NL)t=+oo - {Nr - NL)t=-oo = -2Qw ■ (2) 

The r.h.s. of Eq. ([2]) is the integral over space-time of 
S^jg , where jg represents the anomalous flavor-singlet 
axial current. The physical picture that arises from 
Eq. ([2]) is that in the presence of topological excitations 
such as instantons and sphalerons with a given Qw, a-nd 
starting with a system of quarks with Nr ~ Nl, an un- 
balance between left-handed and right-handed quarks is 
produced. Such an unbalance can lead to observable ef- 
fects to probe topological V- and CV-odd excitations. An 
experimental observable sensitive to local V- and CV- 
odd effects has been proposed in Ref. Recently the 
STAR collaboration presented the conclusive observation 
of charge azimuthal correlations [l^l possibly resulting 
from the CME with local V- and CP-violation. 

The intuitive picture of the CME is as follows. In a 
strong magnetic field B, quarks are polarized along the 
direction of B. Let us suppose that B is along the pos- 
itive z axis (that is conventionally taken as the y axis 
in the context of heavy- ion collisions) . Neglecting quark 
masses, which is a good approximation for u and d quarks 
in the high-T chiral restored phase, the chirality is an 
eigenvalue to label the quarks. Then, right-handed u 
quarks should have both their spin and momentum par- 
allel to B and left-handed u quarks should have their 
spin parallel to B and momentum anti-parallel to B. 
Obviously the same reasoning applies to d quarks. If 
Nr = Nl, then the current that would originate from 
the motion of left-handed quarks is exactly cancelled by 
that of right-handed quarks. If Nr ^ Nl which is ex- 
pected from the anomaly relation on the other hand, 
a finite net current is produced. Therefore, if quarks ex- 
perience a strong magnetic field in a domain where the 
topological transition occurs, a net current is produced 
locally. 

The CME has been investigated in the chiral effective 
model [2^ as well as in the holographic QCD model fl7| . 
The chiral magnetic conductivity is calculated without 
gluon interactions in Ref. [2^ . In (2^ the electric-current 



susceptibility under a homogeneous magnetic field, which 
can be related to the fluctuation of the electric-charge 
asymmetry measured by the STAR collaboration, has 
been computed in the same way. The first lattice-QCD 
study of the CME has been performed by the ITEP lat- 
tice group j28j in the color-SU(2) quench approxima- 
tion. Moreover, the Connecticut group [29| performed 
a lattice-QCD study of the CME with 2-1-1 dynamical 
quark flavors. 

This article is devoted to the study of the two-flavor 
Nambu-Jona-Lasinio model with the Polyakov loop cou- 
pling (PNJL model) in a strong magnetic field. The 
PNJL model has been introduced in Refs. [13, HH to in- 
corporate deconfinement physics into the N JL model [s^j . 
The main addition to the NJL model is a background 
gluon field in the Euclidean temporal direction. The 
background field is related to the expectation value of 
the traced Polyakov loop, which is known to be an or- 
der parameter for the deconfinement transition in a pure 
gauge theory (ssj . There are many theoretical studies re- 
lated to different aspects of the PNJL model; see for ex- 
ample Ref. See Refs. [s^ for a related study within 
the Polyakov-Quark-Meson model, and [s^ for an inves- 
tigation within QCD with imaginary chemical potential. 

We work in the chiral limit throughout the paper, in 
which the definition of the chiral critical temperature has 
no ambiguity. Firstly we focus on the effect of B on chi- 
ral symmetry restoration at finite temperature. As it will 
be clear soon, our results support the role of the external 
B as a catalyzer of dynamical symmetry breaking; the 
critical temperature increases with increasing B. Natu- 
rally the (pseudo)critical temperature for deconfinement 
crossover is less sensitive to the presence of B because 
there is no direct coupling between photons and gluons. 
Hence, the PNJL model predicts that at large enough 
B, a substantial range of temperature will open at which 
quark matter is deconfined but chiral symmetry is still 
dynamically broken. See [13] for a related study. 

Also we shall discuss the effects of a finite chiral chem- 
ical potential /K5 on the phase structure within the PNJL 
model. This /its mimics the topologically induced changes 
in chirality charges TVs = Nr^ — Nl that are naturally ex- 
pected by the QCD anomaly relation. The relevant quan- 
tity in a microscopic picture is rather the total chirality 
charge N^ but for technical reasons it is easier to work in 
the grand-canonical ensemble by treating /ig (see Ref. [ll] 
for an alternative description based on the fiux-tube pic- 
ture) , which is to be interpreted as the time derivative of 
the 6 angle of the strong interactions; fi^ = 0/{2Nf). Be- 
sides the phase diagram from the PNJL model, we com- 
pute quantities that are relevant to the CME, namely, the 
induced electric current density, its susceptibility, and the 
chiral charge density together with its susceptibility. 

This paper is organized as follows; in Sec. HI] we give 
a detailed description of the model we are using. In 
Sees, mil and IIVI we present and discuss our numerical 
results from the model. Finally, in Sec. |V] we draw our 
conclusions. 
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II. MODEL WITH MAGNETIC FIELD AND 
CHIRAL CHEMICAL POTENTIAL 

In this section we analyze the interplay between the 
chiral phase transition and the deconfinement crossover 
at large B using the PNJL model. Here we consider 
two-flavor quark matter in the chiral limit since the chi- 
ral phase transition is a true phase transition only in the 
chiral limit, and then and only then Tc can be identified 
unambiguously by vanishing chiral condensate. The chi- 
ral limit in the two-flavor sector is not far from the phys- 
ical world in which the current quark masses are a few 
MeV, almost negligible as compared to the temperature. 
Moreover, we are interested in studying the situation in 
the presence of chirality charge density. In the grand- 
canonical ensemble we can introduce the chirality charge 
by virtue of the associated chemical potential ^5 in the 
following way. 

The Lagrangian density of the model we consider is 
given by 

C = ij {i^^D^ + Ai57°7') ^ 



G 



{tpijf + {i)i-i^Tijf - u[^, !>, T] , (3) 



where the covariant derivative embeds the quark cou- 
pling to the external magnetic field and to the back- 
ground gluon field as well, as we will see explicitly below. 
We note that /i5 couples to the chiral density operator 

A/'s = '07°7^?/' = VflV'i?^ - V'lV'L, hence ns = ^ 

can develop when /is 7^ 0. The mean- field Lagrangian is 
then given by 

£ = [i-f^D^^ -M + M57°7') V' - T] , (4) 

where M = ~2a with a = G{ij'ip) = G{{uu) + (dd)). 

In Eq. (lU $, $ correspond to the normalized traced 
Polyakov loop and its Hermitean conjugate respectively, 
$ = {1/Nc)TtL and I* = (l/7Vc)TrLt, with the Polyakov 
loop matrix, 



V exp 




(5) 



where ,3 = 1/T. 

The potential term U[^,^,T] in Eq. dH) is built by 
hand in order to reproduce the pure giuonic lattice 
data Among several different potential choices [l^ 
we adopt the following logarithmic form [3l|, [s^ , 



Z^[$,$,T] = T- 



a{T) 



b{T) ln[l - 6$$ + 4($3 + $3) - 3($$)2] I , 



(6) 



with three model parameters (one of four is constrained 



by the Stefan-Boltzmann limit), 



a{T) = ao + ai ( — 



b{T) = h 



0.2 



(7) 



-1.75. 



The standard choice of the parameters reads (3j 

ao = 3.51, ai = -2.47, as = 15.2 , 63 = 

(8) 

The parameter To in Eq. ([6|) sets the deconfinement scale 
in the pure gauge theory, i.e. Tc = 270 MeV. 

We assume a homogeneous magnetic field, JB, along 
the positive z axis. The eigenvalues of the Dirac operator 
can be derived by the Ritus method [i^, which are [l9| : 



[IpI + s/i5Sgn(p,)]' 



(9) 



apart from (the phases of) the Polyakov loop, where s = 
±1, = pi + 2|(7/i?|fc with k a non-negative integer 
labelling the Landau level. 

The thermodynamic al potential in the mean-field 
approximation in the presence of an Abelian chromo- 
magnetic field has been considered in many literatures, 
Ref. [4l| for example. The expression for an electromag- 
netic B can be obtained in the same way. Here we simply 
write the final result; 



G 



f=u,d 



27r 



ask 



dPz r2 



f—it.d s.k 

X lnfl + 3$e-'3"= +3^e-^^' 



dpz 
2tt 



271 



-3/3 



fA^sip) 



(10) 



Here the above definition of ft is different from the stan- 
dard grand potential in thermodynamics by a volume 
factor V. The quasi-particle dispersion ojg is given by 
Eq. ([9]). The spin degeneracy factor is 



ask 




fc = 0, qB>0 , 
fc = 0, qB <0 , 



(11) 



Before going ahead further, one may wonder why we in- 
troduce only one order parameter for the chiral symmetry 
breaking even though the magnetic field breaks isospin 
symmetry. Since Qu ^ Qd, one could suspect that the ef- 
fects of B on (uu) and on (dd) are different. This is not 
the case, however, in the present model in the mean-field 
approximation. As a matter of fact, even in the presence 
of JB 7^ 0, the thermodynamic potential ([TU)) depends 
only on cr cx (uu) + (dd) . This is so only when the four- 
fermion interaction is Eq. ^ with equal mixing of the 
?7(l)yi-symmetric and [/(I) ^-breaking terms. Hence, the 
relevant quantity for the chiral symmetry breaking is just 
one condensate, namely cr, even for B 7^ 0, and there is 



no need to introduce two independent condensates in this 
special case. Even when we consider more general four- 
fermion interaction, the isospin breaking effect is only 
negligibly small. 

The vacuum part of the thermodynamic potential, 
il{T = 0), is ultraviolet divergent. This divergence is 
transmitted to the gap equations. Thus we must spec- 
ify a scheme to regularize this divergence. The choice 
of the regularization scheme is a part of the model def- 
inition and, nevertheless, the physically meaningful re- 
sults should not depend on the regulator eventually. In 
the case with a strong magnetic field the sharp momen- 
tum cutoff suffers from cutoff artifact since the contin- 
uum momentum is replaced by the discrete Landau quan- 
tized one. To avoid cutoff artifact, in this work, we 
use a smooth regularization procedure by introducing a 
form factor /a (p) in the diverging zero-point energy. Our 
choice of /a(p) is as follows; 



A2W 



^2N 



\P 



2N 



(12) 



where we specifically choose N = 10. In the N ^ oo limit 
the above /a(p) is reduced to the sharp cutoff function 
d{A — \p\). Since the thermal part of is not divergent, 
we do not need to introduce a regularization function. 



III. PHASE STRUCTURE WITH CHIRAL 
CHEMICAL POTENTIAL 
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In this section we firstly focus on the system at ^5 = 
and discuss the role of the magnetic field as a catalyzer of 
the dynamical chiral symmetry breaking. We also ana- 
lyze the interplay between chiral symmetry restoration 
and deconfinement crossover as the strength of B in- 
creases. 



A. Results at /is = — chiral symmetry breaking 
and deconfinement 

We analyze, within the PNJL model, the response of 
quark matter to B at = 0. In particular we are inter- 
ested in the interplay between chiral symmetry restora- 
tion and deconfinement crossover in the presence of a 
magnetic field which leads to the so-called chiral mag- 
netic catalysis [ij]- Our model parameter set is 

A = 620 MeV , GA^ 2.2 . (13) 

These parameters correspond to /^r = 92.4 MeV and the 
vacuum chiral condensate {uu)^/^ = —245.7 MeV, and 
the constituent quark mass M = 339 MeV. The critical 
temperature for chiral restoration in the NJL part at B = 
is Tc « 190 MeV. We set the deconfinement scale 
To in the Polyakov loop potential (see Eq. ^) as Tq = 
270 MeV, which is the value of the known deconfinement 
temperature in the pure SU(3) gauge theory. 



FIG. 1. Absolute value of the chiral condensate \ {uu)^^^\ (up- 
per panel) and expectation value of the Polyakov loop (lower 
panel) as a function of T computed at several values of eB 
(in unit of ml). In this model Tc = 228 MeV at fi^ = B = 0. 

In Fig. [l]we plot the absolute value of the chiral con- 
densate (uu)^^'^ (upper panel) and expectation value of 
the Polyakov loop (lower panel) as a function of T com- 
puted at several values of eB (expressed in unit of mj). 
The chiral condensate (uu) and the Polyakov loop $ are 
the solution of the gap equations dVl/da = dVl/d^ = 
in the model at hand. 

Figure [T] is interesting for several reasons. First of all, 
the role of B as a catalyzer of chiral symmetry breaking 
is evident. Indeed, the chiral condensate and thus con- 
stituent quark mass increase in the whole T region as eB 
is raised (for graphical reasons, we have plotted our re- 
sults starting from T = 100 MeV. There is nevertheless 
no significant numerical difference between the T = 
and T = 100 MeV results). This behavior is in the cor- 
rect direction consistent with the well-known magnetic 
catalysis revealed in Ref. [13] and also discussed recently 
in Ref. [4l[ in a slightly different context of the PNJL- 
model study on the response of quark matter to external 
chromomagnetic fields. 

Secondly, wc observe that the deconfinement crossover 
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is only marginally affected by the magnetic field. We 
can identify the deconfinement by the inflection point 
of $ as a function of T. This simple procedure gives 
results nearly in agreement with those obtained by the 
peak position in the Polyakov loop susceptibility, which 
is a common prescription to locate the so-called pseudo- 
critical temperature. Also we can identify the deconfine- 
ment Tc with the temperature at which $ = 0.5. We note 
that Tc in Fig. [T] is the chiral = 228 MeV where the 
chiral condensate vanishes, but not the deconfinement 
in both figures. 

From Fig. [1] we notice that, increasing eB from 4m^ 
to 20m^, the shift of the chiral transition temperature 
AT^ « 20 MeV, while the shift of the deconfinement 
crossover temperature is as small as AT^ k, 5 MeV. 
Hence, the chiral phase transition is more easily in- 
fluenced by the magnetic field than the deconfinement 
as anticipated. Consequently, under a strong magnetic 
field, there opens a substantially wide T-window in which 
quarks are deconfined but chiral symmetry is still spon- 
taneously broken. This result, valid for /Lt5 = 0, does not 
necessarily hold, in general, for /X5 ^ 0, as we will see in 
the next subsection. 



B. Results at /is 7^ — suppression on the chiral 
condensate 

We next turn to the study of the effect of a finite ^jl^ 
on the QCD phase transitions using the PNJL model. 
We recall that /X5 cannot be a true chemical potential 
since its conjugate variable is only approximately con- 
served due to axial anomaly. Nevertheless, because the 
time derivative of the strong 9 angle translates into , as 
explained in the introduction section, itself is a phys- 
ically meaningful quantity. We specifically look into the 
behavior of the critical line for chiral symmetry restora- 
tion, which is well defined in the chiral limit, at differ- 
ing B while keeping /is fixed. This study will be useful, 
among other things, in order to understand the relation 
between the chirality density 715 and /i.s that we compute 
numerically in later discussions. 

One effect of /is 7^ is lowering of the critical tem- 
perature of the chiral phase transition. This is evident 
from the upper panels of Fig. [2] Firstly we discuss the 
case of eB = 5m^ corresponding to the left upper and 
left lower panels of Fig. [21 We see that increasing /^s 
at low T results in slight enhancement of the chiral con- 
densate. As T approaches Tcj however, the chiral phase 
transition at larger /xs takes earlier place below Tc- As a 
result of the coupling between the chiral condensate and 
the Polyakov loop, the deconfinement crossover as shown 
in the lower panels of Fig. [2] is also shifted earlier as //s 
becomes greater. 

In view of the right upper and right lower panels of 
Fig. [5] for large eB ~ lOm^ the /xs-effect on the chiral 
condensate at low T is less visible. This is understood 
from the fact that the chiral magnetic catalysis effect 



is predominant over the minor enhancement due to /is. 
In contrast, as T is increased toward T^, the qualitative 
behavior of the shift in the critical temperature is just 
the same as what we have seen previously for eB = 5to^. 

An interesting prediction from the PNJL model is that, 
at a given value of e_B, there exists a critical /is, above 
which the chiral phase transition becomes first order. 
In the case eB = 5to^ as shown in Fig. [2] the criti- 
cal /is is found between 300 ^ 400 MeV. As a mat- 
ter of fact, at /is = 400 MeV, the chiral condensate 
and the Polyakov loop both exhibit discontinuity at the 
critical temperature. We see that, as compared to the 
/is — 300 MeV case, the slopes of the chiral condensate 
and the Polyakov loop sharply change as a function of T 
for the /is = 400 MeV case. Hence, our data plotted in 
Fig. m suggest the existence of a critical /is in the range 
300 MeV < /x§ < 400 MeV at which the weakly first- 
order transition becomes a true second-order one. The 
phase diagram in the n^-T plane has a tricritical point 
(TCP) accordingly. We will discuss more on the TCP 
in the next subsection. We notice that this picture is 
qualitatively robust regardless of the chosen value of e_B, 
as is already implied from Fig. [2j The first-order phase 
transition in the high-mus and low T region leads to a 
discontinuity in the chirality density as a function of /is . 
This point will be also addressed in some details in the 
next section. 

As a final remark in this subsection we note that, in 
Figs. [T] and [2l the slope of the curve quickly changes at 
some point (at T/Tc = 1 for eB = for example). This 
is because of the second-order phase transition which is 
the case for chiral restoration in the chiral limit. 



C. Phase diagram 

The results we have revealed so far can be summarized 
into the phase diagram in the /is-T plane. In the upper 
panel of Fig. |3]we show the phase diagram at eB = 5toJ. 
In the lower panel for comparison we plot the phase di- 
agram at eB = 15m^. The line represents the chiral 
phase transition. It is of second order for small values 
of /is (shown by a thin line) and becomes of first order 
at large /is (shown by a thick line). The location of the 
TCP on the phase diagram depends only slightly on eB, 
while the topology of the phase diagram is not sensitive 
to the magnetic field. 

The general effect of /^s is to lower the chiral transition 
temperature. One may argue that the critical line can 
hit T = eventually at very large /xs, though the PNJL 
model is of no use at such large /xs because the ultraviolet 
cutoff causes unphysical artifacts. The locations of the 
TCP are estimated from the PNJL model as 

(/is, T) w (400 MeV, 200 MeV) , for eB = 5ml , (14) 
(/is,r) « (370 MeV, 200 MeV) , for eB = 15m^ . (15) 
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FIG. 2. Absolute value of the chiral condensate |(itit)^''''| (upper panel) and expectation value of the Polyakov loop (lower 
panel) as a function of T computed at several values of eB (in unit of mj) and /15 (in unit of MeV). 



IV. CHIRALITY CHARGE, ELECTRIC 
CURRENT, AND SUSCEPTIBILITIES 

In this section we show our results for quantities rele- 
vant to the Chiral Magnetic Effect (CME). We numeri- 
cally compute the chiral density 715 and its susceptibility 
X5 as a function of and eB. Also we calculate the cur- 
rent density j'3 along the direction of B and its suscep- 
tibility xj. Finally, we use the result n5(/^5) to evaluate 
jg as a function of 71,5 . 



A. Chirality density and its susceptibility 

The axial anomaly relates the topological charge Qw 
to the chirality charge A^s with TVs = 't-s^ where V = 
L^LyLz is the volume of topological domains. We can 
read from 

n5 = . (16) 



It is useful information to relate 715 and /is for various 
temperatures and magnetic field strength. In the next 
section we will use the results of ?T.5(/i5) to express the 
current density as a function of the chirality density. 

The relation between 715 and /is can be found analyti- 
cally only in simple limiting cases [l9j . In general one has 
to determine it numerically using an effective model. We 
show ns(/i5) for eB = 5m^ at three temperatures around 
Tc in Fig. 21 The qualitative picture is hardly modified 
even if we change the magnetic field. 

From Fig. [3] we can read the critical temperature at 
/is = that is Tc = 228 MeV. At temperatures well be- 
low Tc, as seen in the T = 160 MeV case in Fig. H the 
discontinuity associated with the first-order phase transi- 
tion with respect to (uu)^^^ and $ is conveyed to the rela- 
tion ns(/is)j which is a typical manifestation of the mixed 
phase at critical /is. Naturally, as T gets larger, the chi- 
rality density as a function of /is becomes smoother, since 
the chiral phase transition is of second order at higher T 
as is clear from Fig. [S] 

It is interesting to compute the chirality charge sus- 
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FIG. 3. (Upper panel) Phase diagram in the fi5-T plane ob- 
tained at eB = 5m^. The thin line represents a second-order 
chiral phase transition and the thick one a first-order transi- 
tion. Below the line, chiral symmetry is spontaneously bro- 
ken, while chiral symmetry is restored above the line. The la- 
bel "TCP" denotes the tricritical point. (Lower panel) Phase 
diagram for eB = 15m^. 



[fm-3] 
10 



8 



in 
C 





/' ^ 


r= 160 MeV 


/ / 


r=220MeV 


// 
/ / 


r= 240 MeV 






/ 




^' / 1 
. / / 
y j 







100 200 300 400 500 600 
Us [MeV] 



FIG. 4. Chirality density 715 (in unit of fm ■^) as a function 
of /is (in unit of MeV) at eB — 5m^ for several values of T. 



ceptibility X5, as well as 715, defined as 

1 d^n 



X5 



13V d^il 



(17) 



where (3 = 1/T and V the volume. We note that this 
definition of the susceptibility is different from that in 
Ref. [l^l by It should be mentioned that we take 
a numerical derivative to compute X5 including implicit 
dependence in $ and cr. In Fig.|5]we plot ^ ^ function 
of eB for several /^s values. The upper panel corresponds 
to /is = 0, middle one to ^(5 = 200 MeV, and lower 
one to fjL^ = 400 MeV. For completeness, in the right 
panels of the same figure [5] we plot the chiral condensate 



\1/3| 



for the same T and same ^(5. The oscillations in 



X5 are artificial results because of the momentum cutoff 
A. As shown in Ref. [4lj, choosing a regulator which 
is smoother than used in this work, the oscillations of 
the various quantities could be erased. The qualitative 
picture is, nevertheless, unchanged even with a different 
regulator. For this reason we do not perform a systematic 
study here on the cutoff scheme dependence. 

A notable aspect is the suppression of the chirality- 
charge fluctuations at large T and large eB. This is 
evident, for example, in the result with /is = and 
T = l.lTc in Fig. [S] As long as eB is small, xs is 
a monotonously increasing function of eB as expected 
naively. When eB reaches a critical value around 20m J, 
however, xs has a pronounced peak and then decreases 
with increasing eB, which is a result of mixture with di- 
verging chiral susceptibility at the chiral phase transition. 
It should be mentioned that xs at /is = (as shown in 
the upper left panel of Fig. [5]) does not diverge at the 
critical eB since the mixing with the chiral susceptibility 
is vanishing due to /is = 0. This behavior below and 
above the chiral critical point can be easily understood 
in terms of the chiral symmetry breaking by virtue of 
the magnetic field. As a matter of fact, ai T > Tc the 
chiral condensate stays zero identically as long as eB is 
small enough, leading to zero quasiparticle masses. Once 
eB exceeds a critical value, the chiral symmetry is bro- 
ken spontaneously even at high T (see the upper right 
panel of Fig. [5|) and the quasiparticle masses can then 
jump to a substantially large number then. Such dy- 
namical quark masses result in appreciable suppression 
of the chirality-charge fluctuations. As it will be shown 
in the next section, this interesting and intuitively un- 
derstandable effect appears in the current susceptibility 
as well. 



B. Current density and its susceptibility 

The current density J3 (and its susceptibility as well) is 
the most important quantity to compute for the Chiral 
Magnetic Effect [l^l- It corresponds to the charge per 
unit volume that moves in the direction of the applied 
magnetic field in a domain where an instanton/sphaleron 
transition takes place, which causes chirality change of 
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quarks. The current has been computed analytically in 
Ref. [1^ in four different ways. 

To compute the current density along the magnetic 
field, i.e. = q{ip"f3ip), we follow the common proce- 
dure to add an external homogeneous vector potential 
along the magnetic field, A3, coupled to the fcrmion field. 



Then, 



J3 



(18) 



^3=0 



The derivative of the thermodynamic potential in the 
presence of a background field is computed in the follow- 
ing way. The coupling of quarks to A3 is achieved by 



9 



shifting pz in Eq. ((TO)) as Pz ^ Pz + qfA^. After putting 
regularization in the momentum integral with an uhra- 
violet cutoff A (we know that the current is uhraviolet 
finite, hence the choice of the regularization method docs 
not affect the final result) we change the order of the mo- 
mentum integral and the derivative with respect to A3. 
Then we make use of the following replacement, 



d 
dAs 



If 



to obtain, 



J3 



f=u,d 



llfBl 
27r 



s.k 



dpz 



^ dpz d 

— — [wsb) 



(19) 



2tt dpz 



(20) 

The ellipsis represents irrelevant matter terms. After 
summing over the spin s, the contribution of the inte- 
grand from the Landau levels with n > turns out to 
be an odd function of pz- Therefore, only the lowest 
Landau level gives a non-vanishing contribution to the 
current and we get from the surface contribution [l9j . 



j3 = iVc E 

f=u,d 



27r2 



67r2 



(21) 



which is certainly ultraviolet finite as it should be. Gen- 
erally speaking we should utilize a gauge-invariant regu- 
larization. Nevertheless the above ([21]) indicates that a 
naive momentum cutoff can reproduce a correct expres- 
sion for the anomalous chiral magnetic current. 

The current density as given by Eq. (|2ip does not de- 
pend on quark mass explicitly, and on temperature ei- 
ther. The reason is that the current is generated by the 
axial anomaly and it receives contributions only from the 
ultraviolet momentum regions (as the above derivation 
shows), and so it is insensitive to any infrared energy 
scales. Also, the Polyakov loop does not appear explic- 
itly in Eq. (PT|) . This is easy to understand; the Polyakov 
loop is a thermal coupling between quark excitations and 
the gluonic medium, and thus the Polyakov loop only en- 
ters the thermal part of fl. Since the current originates 
from the anomaly, however, the thermal part of just 
drops off for the current generation. The effect of the 
Polyakov loop will appear implicitly through the relation 
between ^5 to 715. 

To confirm that our numerical prescription works well, 
we have computed J3 by means of Eq. ([TSl) with given in 
Eq. (jlOp . In Fig.[6]we show the results from our numerical 
computation as a function of ^5. In the figure we have 
plotted the normalized current. 



J3 



5^0 e^B 
67r2 



J3 



(22) 



with a choice oi fiQ ~ 1 GeV, which we defined so to make 
the comparison transparent at a glance. In Fig. [5] the 
dashed line represents at T = 160 McV; on the other 
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FIG. 6. (Upper panel) Normalized current density, js = 
(5/io e^_B/67r^)~^j3 with fj,o — 1 GeV, as a function of 
at two different temperatures (below and above Tc). (Lower 
panel) Current density as a function of ns for eB — brrii 
computed at three different temperatures. 



hand, the dot-dashed line the case at T = 240 MeV. We 
notice that our numerical results are perfectly in agree- 
ment with Eq. pT|) . We conclude that our numerical 
procedure correctly reproduces the expected dependence 
of on ^5 with the correct coefficient insensitive to in- 
frared scales regardless of whether T is below or above 

The result shown in the upper panel of Fig. [5] gives 
us confidence in our numerical procedure but the figure 
itself is not yet more informative than Eq. ((2T|) . We ex- 
press now j3 as a function of using Eq. (|2ip and the 
results discussed in the previous section. The result of 
this computation is shown in the lower panel of Fig. [HI in 
which we plot the (not normalized) current density (mea- 
sured in fm"'^) as a function of 715 (measured in fm"''), 
at eB — 5to2 and at three different temperatures. 

From Fig. [6] we notice that, at a fixed value of ns, the 
larger the temperature is, the smaller becomes. This 
seemingly counter intuitive result is easy to understand. 
As a matter of fact, as the temperature gets larger, the 
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FIG. 8. Subtracted current susceptibility with a smoother 
regulator with TV = 5 in Eq. (fT2)) . 



This is in contradiction with the gauge invariance, which 
requires the above susceptibility to vanish because the 
current susceptibility is nothing but the 33-component 
of the photon polarization tensor at zero momentum. 
Therefore, in order to deal with the physically meaning- 
ful quantity, we subtract the vacuum part from the above 
equation and compute, 



> 

CD 
O 

> 

=0, 




XJ = Xj(M5,5,T) - xj(/i5,0,0) 



(24) 



eB/m/ 



FIG. 7. Subtracted current susceptibility, PVxj, as a function 
of eB (in unit of m^) for several different values of T (in unit 
of Tc = 228 MeV) and /is (measured in MeV). 



corresponding fi^ for a given should decrease because 
of more abundant thermal particles at higher tempera- 
ture. Since js depends solely on /ts, a higher temperature 
requires a larger 715 to give the same J3. 

Besides J3, another interesting quantity is the current 
susceptibility defined by 



XJ 



01) 



{.13 



1 d^n 



PV dAl 



(23) 



If we naively use the above definition ([23| for the cut- 
off model like the PNJL model, xj is non-zero propor- 
tional to even at T = i? = as discussed in Ref. 271 



to fulfill the requirement that photons are unscreened at 
T = B = regardless of any value of 

We plot our results for (3Vxj as a function of eB in 
Fig. [7] at /is = (upper panel), /X5 = 200 MeV (middle 
panel), and /X5 = 400 MeV (lower panel). The oscilla- 
tions in the susceptibility behavior are an artifact of the 
momentum cutoff. In these plots Tc — 228 MeV denotes 
the critical temperature for chiral symmetry restoration 
at /t5 = i? = 0. 

Let us first focus on the case at /X5 = 0. At T = 
and T = 0.95Tc the system is in the broken phase with 
(uu) 7^ over the whole range of eB. On the other 
hand, at the temperature T = l.lTc, the system is in 
the chiral symmetric phase for eB smaller than a critical 
value. There is a phase transition from the symmetric to 
the broken phase with increasing eB. This transition is 
driven by the presence of the magnetic field as the catal- 
ysis of chiral symmetry breaking, as mentioned before. 
The effect of the phase transition leads to a cusp in the 
susceptibility x,/ as a function of eB. We also notice that 
there seems to exist a range in eB in which xj < 0. This 
is a mere artifact of the momentum cutoff, which causes 
unphysical fluctuations in x.i- The qualitative picture is 
similar also at /i5 7^ 0. 

To distinguish physically meaningful information from 
cutoff artifacts, we have computed XJ using a smoother 
regulator with iV = 5 in Eq. (fT2|) . We have readjusted 
the NJL parameters to keep the physical quantities (/^r 
and (uu)) unchanged. Figure [S] is the result in which 
oscillations are suppressed and xj > for any T and B. 
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In view of Figs. [7] and [8] we can conclude that it is 
certainly a physical effect that the chiral phase transi- 
tion critically affects the susceptibility XJ well as Xb- 
As shown in Ref. [l^] the susceptibility difference be- 
tween the longitudinal and transverse directions has an 
origin in the axial anomaly and is insensitive to the in- 
frared information. Nevertheless, xj (and transverse Xj 
too) should be largely enhanced near the chiral phase 
transition through mixing with the divergent chiral sus- 
ceptibility, which is not constrained by anomaly. Such 
enhancement in xJ would ease the confirmation of the 
CME signals at experiment. 

V. CONCLUSIONS 

In this article we have considered several aspects re- 
lated to the response of quark matter to an external 
magnetic field. Quark matter has been modelled by the 
Polyakov extended version of the Nambu~Jona-Lasinio 
(PNJL) model, in which the QCD interaction among 
quarks is replaced by effective four-fermion interactions. 
In the PNJL model, besides the quark-antiquark con- 
densate which is responsible for the dynamical chiral 
symmetry breaking in the QCD vacuum, it is possible 
to compute the expectation value of the Polyakov loop, 
which is a relevant indicator for the quark dcconfinemcnt 
crossover. 

In our study, we have firstly focused on the effect of a 
strong magnetic field on chiral symmetry restoration at 
finite temperature. Our results show the effect of the ex- 
ternal field as a catalyzer of dynamical symmetry break- 
ing. Moreover, the critical temperature increases as the 
strength of B is increased. This behavior is in agree- 



ment with the previous studies on magnetic catalysis in 
NJL-like models. 

We have also discussed the effects of a chiral chemical 
potential, /i5, on the phase structure of the model. The 
chiral chemical potential mimics the chirality induced by 
topological excitations according to the QCD anomaly 
relation. Instead of working at fixed chirality iVs, we 
have worked in the grand-canonical ensemble introducing 
/i5, i.e. the chemical potential conjugate to N^,. Besides 
the phase diagram of the model, summarized in Fig. |3l 
we have computed several quantities that are relevant for 
the Chiral Magnetic Effect (CME). That is, we have com- 
puted the current density j3 and its susceptibility xJ as 
well as the chiral charge density ris and its susceptibility 

As a future project it is indispensable to extend our 
analysis to the 2 -f 1 fiavors and tunc the PNJL model 
parameters to reproduce the correct Tc and thermo- 
dynamic properties, which would enable us to make 
a serious comparison with the dynamical lattice-QCD 
data d^l, and furthermore, it would be possible to give 
more pertinent prediction on the physical observables. 
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